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ABSTRACT 

The nature of solar wind (SW) turbulence below the proton gyroscale is a topic that is being inves- 
tigated extensively nowadays, both theoretically and observationally. Although recent observations 
gave evidence of the dominance of Kinetic Alfven Waves (KAW) at sub-ion scales with uj < uJd, other 
studies suggest that the KAW mode cannot carry the turbulence cascade down to electron scales and 
that the whistler mode (i.e., u > uJd) is more relevant. Here, we study key properties of the short 
wavelength plasma modes under limited, but realistic, SW conditions, typically /3i > /3e ~ 1 and for 
high oblique angles of propagation 80° < OkB < 90° as observed from the Cluster spacecraft data. The 
linear properties of the plasma modes under these conditions are poorly known, which contrasts with 
the well-documented cold plasma limit and/or moderate oblique angles of propagation (OkB < 80°). 
Based on linear solutions of the Vlasov kinetic theory, we discuss the relevance of each plasma mode 
(fast, Bernstein, KAW, whistler) in carrying the energy cascade down to electron scales. We show, 
in particular, that the shear Alfven mode (known in the MHD limit) extends at scales kpi > 1 to 
frequencies either larger or smaller than iOd, depending on the anisotropy fcy/Zcj^. This extension 
into small scales is more readily called whistler (w > ojd) or KAW {oj < uJd), although the mode is 
essentially the same. This contrasts with the well-accepted idea that the whistler branch develops 
always as a continuation at high frequencies of the fast magnetosonic mode. We show, furthermore, 
that the whistler branch is more damped than the KAW one, which makes the latter the more relevant 
candidate to carry the energy cascade down to electron scales. We discuss how these new findings 
may facilitate resolution of the controversy concerning the nature of the small scale turbulence, and 
we discuss the implications for present and future spacecraft wave measurements in the SW. 
Subject headings: turbulence, solar wind, heating, whistler, KAW 



1. INTRODUCTION 

The nature of the SW turbulence below the ion scale 
(typically pi ~ 100 km, corresponding to an observed 
frequency in the spacecraft frame of fsc ^ 0.5 Hz) has at- 
tracted considerable interest in the space and astrophys- 
ical communities in recent years. This has been encour- 
aged in particular by recent observations from the Clus- 
ter mission that provided the most complete and detailed 
picture of the SW turbulence cascade from magnetohy- 
drodynamic (MH D) scales (L » pj) to electron scales 
(L ~ Pe ~ 1 km) (iS ahraoui et al.|20"09 iKiyani et al.|2009 ' 



Alexandrova et arf 2009: C hen et al. ''2011;~Sahr aoui et| 
al.| |2010ap . Pet ermining the nature and properties (e.g., 
scaimg, anisotropy) of the turbulence at small scales is 
indeed a crucial point to understanding the problems of 
energy dissipation and heating, particle acceleration, and 
magn etic r econnec tion in space and astrophysical plas- 



mas ( Schekochihin et al. 20091. Recent Cluster obser 



vations provided clear evidence that SW turbulence cas- 
cades below the ion scale pi down to the electron scale 
pe where dissipation becomes importan t and the spec- 
tra steepen to ^ /sc~", with a > i (Sahraoui et al. 
2009 2010a I. The spectrum thus formed was termed the 
dissipation range while the range of scales between pi 
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and Pe has been termed the dispersive range in reference 
to the dispersive nature of the plasma modes at those 
scales [Stawicki et al.„2001j ). In a case study, Sahraoui 



et al. (2010a) pertorrnied a detailed analysis ot the en 
ergy cascade from MHD to sub-ion scales using a multi- 
point measurement technique called fc-filtering ( Pingon & 
'LefeuvreT991 : Sahraoui e{^al.|2003al|2010b||Narita et aT 
|2010a b, Tjulin et al. 2005| ). The results showed clearly 
that: i) I'he magnetic turbulence is strongly anisotropic 
{kj_ » fc||) down to the observed scale kj_pi ^ 2; ii) The 
ca scade is consistent wi th K AW turbulence as proposed 
in [Howes et al] ( l2008a[ ) and [Schekochihin et al.| pOOQ l 
with frequencies m the plasma rest frame ui O.lwci al- 
though the frequencies in the spacecraft frame reached 
20a;ci; iii) The turbulence undergoes a transition range 
near the ion scale pi characterized by a steepening of 
the spectrum from kj^'^ to fcj"*'^, which has been in- 
terpreted as due to Landau damping of magnetic energy 
into ion heating ^ Howes et al.|2008b I . The remaining en- 
ergy cascades following a power-law ~ fTc^'^ down to the 
electron scale where the energy was suggested to dissi- 
pate into electron heating. The wave number spectra at 
those frequencies could not have been directly measured 
using the fc-filtering technique due to the limitation im- 
posed b y the satellites separati ons, which are larger than 
100 km ( [Sahraoui et"ar[[2010c[ ). However, ID wavenum- 
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ber spectra can be inferred indirectly by means of the 
Taylo r frozen-in-flow approximation provided that it is 
valid ([Bale et al. [2005 Sahraoui et al.|2d09 Alexandrova 



et al. |2UU9p . Sirnilar analyses of magnetic turbulence at 



MMD and sub-ion scales using the fc-filtering technique 
have been carried out recently and have confirmed most 



of these results (|Narita et al.||2011 


Sahraoui et al.||2011 


Narita et al.| 


2010a|b||^ 



erally polarized between those who advocate for the rel- 
evance of Kinetic Alfven Wave (KAW) turbulenc e to ex- 



plain observations of small scale SW turbulence ( Howes 
erar] [2008a|bl |Schekochihin et al.| [20091 [Howes et aT 
2011ap and those who believe in the necessity that 
whistler (or anoth er type) o f turbulence should exist at 
small scal es (Stawicki et al.,,2001 ; Gary & Smith 2009; 



Podesta et aTr' 2010p . indeed, [Podesta et al 



have used a toy model of energy cascade where the lin 
ear damping rate of the KAW waves is considered and 
showed that KAW cascade cannot reach the electron 
scale, because the energy flux of the cascade vanishes at 
scales kpi ^ 20. However, recent GyroKinetic (GK) sim- 
ulations, which self-consistently contain kinetic damping 
of the plasma modes, did not confirm that conclusion and 
showed rather that the KAW mod e can carry turbulenc e 



cascade down to electron scales (Howes et al. 2011a 



Recent 2D PIC simulations showed also the formation of 
a powe r law spectrum kj^'^ below t he electron inertial 



length ( Camporeale fc Burgess||2011 1. In th e framework 
of incompress ible Electron-MHD turbulence [Mayrand fc| 



Galtier (2010) showed that the magnetic energy spectra 
should follow a power law k~^^/'^ at scales smaller than 
de (the electron inertial length), which they proposed to 
explain the steep ening of the spectra near the electron 



scale reported in Sahraoui et al.| (2009 2010a). Note 
however that this nuid model is non-dissipative and does 
not consider any damping of the turbulence via kinetic 
effects, which are important in the dispersive and the 
dissipation ranges. 

Here we derive and analyze the linear solutions of the 
Vlasov-Maxwell equations to shed light on new aspects 
of small scale kinetic plasma modes. Based on the ob- 
tained damping rates of the waves and their magnetic 
compressibilities we discuss the relevance of each mode 
to carry the energy cascade in the SW down to electron 
scales. Unlike the previous papers that had dealt with 
this problem, in this article we focus only on the prop- 
erties of the plasma modes under a restricted range of 
SW parameters, namely hot plasma with /3i > /3e ~ 1 
and high oblique angles of propagation, typically 80° < 



nique (Sahraoui et al.|2010a 


Naritaet al.|2011| 


Sahraoui 


et al.||2011 


). I'his range ot parameters has not been ad- 



solving part of the ongoing controversies on the problem 
of energy cascade and dissipation in the dispersive range. 
Before studying in detail the plasma kinetic modes, in 



1 [Narita et 
plasma modes 



al.| ( |2011^ have claimed to have observed several 
other than the KAW, based on the random spread 
of the observed dispersion relations. However, given the absence of 
any statistical error bars in those observations, one cannot reliably 
interpret those results. A detailed discussion of these important 
issues will be given elsewhere. 



the next section we will first derive the linear solutions 
of the two-fluid theory in the cold and hot plasma cases. 
Although one may question the validity of the fluid de- 
scription of hot plasmas (we do not address that issue 



here, see |Krauss-varban et al.| (|1994P; |Howes | (12009])), 

loluti 



we believe nevertheless that deriving the fluid solutions 
is a necessary step to better understand the more com- 
plex solutions given by the kinetic theory. This com- 
plexity is essentially due to the damping of the modes 
by kinetic effects, which makes them hard to trace at all 
scales, in addition to the rising of several mode conver- 
sions ^Krauss-varban et al.|l994| |Li fc Habbal|200l| . As 

The 



we will show, the reduced-two-fluid theory indeed helps 
tracking the various plasma modes and understanding 
the important changes in their properties depending on 
the plasma /3 and the obliquity of the waves. It allows us 
also to derive analytically various asymptotic properties 
of the waves that will be used to discuss new properties of 
the Alfven and the whistler modes. In the last section we 
discuss the results and their implications on spacecraft 
measurements of waves in the SW. 

2. LINEAR SOLUTIONS OF THE REDUCED-TWO-FLUID 
MODEL 



We 



use the re duced- t wo-fluid theory as described 
(2003b I, whose equations are the fol- 



in Sahraoui et al 
lowing 

dtHs + V.(nsVs) = 



VP, 



P. 



{msnsY- 



= 



dtB = -V X E 
V X B ^0 J 



(1) 

-(E + v, xB) (2) 

(3) 

(4) 
(5) 

where the subscript s denotes ions and electrons, 7, is 
the polytropic index and J = qsUs'Vs is the elec- 
tric current. The system of equations (1 5) is essentially 
equivalent to the classical two-fluid moaclwhcn the non- 
relativistic and quasi-neutrality assumptions are consid- 
ered. This translates into neglecting both the longitudi- 
nal and the transverse components of the displacement 
current, which yields the simplified form of equation [Sj 
This system has the advantage of ruling out the three 
high frequency modes of the two-fluid theory: the optic 
modes (with ~ c, the speed of light) and the Langmuir 
mode uj ojpe. However, it still includes the electron in- 
ertia, which allows us to describe properly the whistler 
mode. The system has therefore three eigenmodes whose 
dispersion relations can be derived analytically. The 
main simplification with respect to the complete kinetic 
system that will be studied in the next section concerns 
the introduction of polytropic closure equations that ex- 
cludes any resonant phenomenon, in particular any effect 
of wave damping. 

To derive the linear solutions we linearize equations ([T]- 
pi) by assuming that Ug — uq + Sug , Pg = Pq + SPg , B = 
Bq + S'B. If, furthermore, we assume that the pertur- 
bations vary as e"-''-'^*"*' '"^ one can obtain (after some 
calculations) the following set of equations 



(5E = BoM, ■ 6v, 



(6) 
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5E - BoM, ■ 5^, 

So 



where 



■ ^ f 

3 1- 



-Jli 



UJLOr 



Jle 



UJU), 
^the 



ekSk + G 

i 

GkBk + G 



(7) 
(8) 

(9) 
(10) 
(11) 



where I is the identity tensor, G = GyGx — GxGy and ek 
is the unit wave vector assumed to be in the XZ plane 



Gk = k/fc = 



Bo = BqBz, uIcs — eB/ms, Vths and Vas are the back- 
ground magnetic field, the gyropulsation, the thermal 
and the Alfven speeds of the particle s, respectively. Re- 
grouping equations (9pT| we obtain the final equation 




M, ^ - ^ - M 



• = 



which admits solutions only if 



Det 



Me - M, - Me • M • M, 



= 



(12) 



Equation [12] yields the final dispersion equation of the 
reduced-two-fluid model 



(1 + /ie.)(l + 2/3,7)^' + (1 + /4)7^' 



-f ii,,f{l + c2 + pa)X^ + c'(l + Me^) 



+ /le»(l+Me»)(l + 2/3,7) 



y4 



(1+Me») 



/5^'(1 + Me,) + W^^^e^{l + A'ez)^' + ^^l^X^ 



0(13) 



where X = kpi, Y = cj/wci, Pi = Vthjood^ Pi = 
Vthi^ /Va^, c = cos^kB, A*ei = rnf,/mi and 7 = 7e -I- 7i is 
the total polytropic index. The solutions of equations [13] 
are plotted in Figs. [T][2] in low and high Pi plasmas for 
the given angles of propagation and SW parameters. 

At low frequencies (w < Wci) we found that the slow 
magnetosonic mode in both low and high /? has an 
asymptotic frequency 

LOci cos 6'kB 

At high frequencies (w > ooci) both hot and cold plasmas 
have two modes with different asymptotes. Fig.]2]shows a 
similar dispersion curve in log-log scale. One can see that 
the first mode is connected at low frequency (w < Ud) to 
the shear Alfven mode (known in MHD) and at kdf. >> 1 
develops an asymptote uj_ given by 

W_ ^ CJce cos 0kB 

It is worth recalling here that the asymptote a;_ is 
usually attributed to the whistler mode, which in this 
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Fig. 1. — Dispersion relations of the reduced-two-fluid model in 
the hot (red) and cold (black) plasmas with = Te, 7; = 5/3 
and 7e = 1. The vertical dashed lines correspond to the electron 
inertial length kde = 1 for the two values of Pi. The insert in the 
top panel is a zoom on the very small scales of the plot to show 
the asymptotes ujj^ and u)— . 
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Fig. 2. — Dispersion relations on a log-log scale of the reduced- 
two-fluid model in the hot (red) and cold (black) plasmas showing 
the connection between high and low frequency modes (the same 
description as in Fig. [l] applies). 

case is connected at low frequency to the Alfven mode 
and not to the fast magnetosonic mode. Therefore, in 
the following we will refer to this mode as the Alfven- 
whistler mode for the sake of clarity. We recall here that 
the designation of slow and fast modes is based upon 
the ordering in the phase speeds of the modes. While 
this ordering remains valid in fluid theories at low and 
high frequency, it can be totally violated in kinetic theory 



(|Krauss-varban et al. 1994 


Li & Habbal 


snows also that the second 


Ub' mode is 



connected at low frequency (both in high and low 
to the fast magnetosonic mode. Therefore, we will refer 
to this mode as the classical-whistler mode. In the low 
Pi limit the classical-whistler mode undergoes significant 
curvature near uj^ at all angles of propagation before 
converging toward the magnetosonic dispersion cj kCg , 
where Cs = y^T^j^e + liVthi is the sound speed and Tgi = 
Te/Ti. From equation 13 one can easily show that in the 
cold case = 0) the Asymptote" aj+ is given by 



L0+ ^ \/ CJce^ COS^ 6'kB + WciWce 

which has two interesting limits (for cjci << i-^ce) 



lim bjj. 



lim 



lim W- = 



The first limit shows that the branch tends to- 
ward the known whistler asymptote uj+ ~ uJce in quasi- 
parallel propagation as shown in Fig. [T] (top panel), while 
the second limit shows that aj+ tends toward the lower- 
hybrid frequency lulh ~ \/^ci^ce. ~ 42LiJci for quasi- 
perpendicular propagation. In contrast, and as can be 
seen in Fig. [T] the asymptote of the Alfven-whistler 
mode continues decreasing as 6'kB — ^ 7r/2 in both low and 
high /3i and becomes even smaller than w^i at all scales 
for the angles of propagation 



For the real mass ratio /iei=l/1836 one obtains Ocrit. = 
89.97°. This implies that at quasi-perpendicular propa- 
gation (i.e., 6'kB > Ocrit.) there are two asymptotes below 



uJci'. one at extremely low frequency, which is the slow 
magnetosonic mode discussed above (wq '■^ ojci cos (^kB) 
and the other one is the Alfven-whistler mode with the 
asymptote a;_ ^ ujce cos 0kB ■ This is the first important 
conclusion of this part. It is worth recalling that these 
two modes were found to be linear solutions of the Hall- 



MHD equa t ions in the in compressible limit (Sahraoui 



et al. 



2007 



Galtier 



2006 ), and become degenerate at 

large (iVLHL)) scales. 'I'hey were referred to respectively 
as the Alfven and the whistler modes. In such incom- 
pressible fluid theories, the classical-whistler mode is in- 
deed ruled out by the incompressibility assumption (i.e., 
/3 — >■ oo). All turbulence theories built in the framework 
of those incompressible fluid models are thus bas ed on the 



slow and Alfven-whistler modes discussed here ( Galtkr 
,2006, 2008). The second important finding from FigTfflis 



the major change that the classical-whistler mode under- 
goes in high Pf. when approaching quasi-perpendicular 
propagation the curvature near tj+ disappears and the 
dispersion curve tends toward the magnetosonic mode 
io ^ kCs- As we will show in the next section, when 
kinetic effects are considered, the classical-whistler mode 
does not simply extend above ujd , but rather splits there 
into different ion Bernstein modes. Becasue the slow 
magnetosonic mode is strongly damped in high /3 plasma 
by kinetic effects, the only relevant mode to carry the en- 
ergy cascade down to the electron scale is the new Alfven- 
whistler mode. 

3. LINEAR SOLUTIONS OF THE VLASOV-MAXWELL 
EQUATIONS FOR HIGH OBLIQUE PROPAGATION 

In this section we solve numerically the line ar Maxwell- 



Vlasov equations using the WHAMP code (Ronnmark 
19821, and compare to the previous linear solutions of 



the reduced hot two-fluid theory. We assume Maxwellian 
distributions functions of electrons and ions and use re- 



alistic SW parameters reported in Sahraoui et al. (2009). 
Part of the obtained dispersion relations and damping 
rates are given in Fig. (the slow magnetosonic mode 
was found to be heavily damped and thus not plotted 
here, nor is that mode studied in the rest of this pa- 
per). The plot shows that the fast magnetosonic modes 
at high oblique angles undergo resonances at the ion cy- 
clotron frequency, which contrasts with the more known 
fast mode in quasi-parallel propagation or low /^^.This 
result is conflrmed by the analysis of the wave polariza- 
tion. In Fig. |4] we plot the phase of the electric field 
component Ey whose sign gives the sense of polariza- 
tion: Arg{Ey) > (resp. < 0) for right (resp. left) hand 
polarized waves. We see a clear transition from right 
to left hand polarization of the fast mode near kpi 1 
(corresponding to a; ^ Wci)- 

In Fig. [3] we find also the kinetic counterpart of the 
Alfven-whistler mode derived from the hot-two-fluid the- 
ory. The Alfven-whistler mode extends the shear Alfven 
mode to small scales as can be seen in the ins ert of Fig. [3 
(th i s resu lt was also found in Fig. 1 of Krauss-varban et 
al. (1994) although not discussed in detail in that pa- 
per)^ The mode becomes dispersive at scales kpi > 1 
and develops frequencies larger (resp. smaller) than LOd 
for 6'kB < 88° (resp. > 88°) up to the scale kp, < 10. We 
refer to the branches oj < Ud and uj > Wci respectively 
as the KAW and the whistler modes. Note that the limit 
UJ = UJci is reached at different spatial scales depend- 
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Fig. 3. — Linear solutions of the Maxwell- Vlasov equations: dis- 
persion relations (blue) and damping rates (red) for the angles of 
propagation 80° < 6kB 1^ 89°. The insert is a log- log plot of the 
same dispersion relations to show the connection between low and 
high frequency modes. 

ing on the value of the angle 0kB • Figure [3] shows that 
in contrast to the fast magnetosonic mode, the Alfven- 
whistler mode has a right hand polarization at all scales 
and does not undergo significant change near kpi ^ 1. 
Nevertheless, the Alfven-whistler modes show features 
related to wave-particle resonances near the harmonics of 
ions (because they are not strictly circularly polarized). 
This can be seen in Fig. [Sl which shows the enhancement 
of the damping of the Slfven-whistler modes near the 
harmonics lon = Nujd — fcyVt/^ (a similar observation 
can be made on Fig. M at other angles of propagation) . 

Fig. [6] shows the Slfven-whistler solutions extended 
to high frequencies and small scales. We observe that 
the damping of the Alfven-whistler mode becomes more 
important when departing from 0kB ^ 90° toward less 
oblique angles. For = 89.99° the solution extends 
down to the electron gyroscale Pe where the damping rate 
remains smalj^ 7/0;^^ 0.4. For less oblique angles the 
Alfven-whistler mode develops frequencies higher than 
uJci but they are subject to stronger damping. This 
can be seen clearly in Fig. [7) which shows the damp- 
ing rate in one period of each wave mode. Figure [7] 
(top panel) shows that the Alfven-whistler modes are 
abruptly damped at kpi ~ 1. At smaller scales the most 
oblique modes are the least damped. The linear damping 
rate may thus play a role of a "filter that lets pass" only 
very oblique modes at small scales. This may explain 
the high oblique modes frequently observed in the SW 
(Sahraoui et al.|2010a||Narita et al. 120111 ISahraoui et al. 

MT|. ' " ' " 

A~similar extension of the fast mode solutions to fre- 
quencies higher than ojci is shown in Fig. |8] We observe 



This solution has been used in [Sahraoui et al.| | |2009| to inter- 
pret SW observations. 
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Fig. 4. — Phase of the electric field component Ey of the Alfven- 
whistler and the fast modes for the angles of propagation 80° < 
^'kB 5: 89°: it is positive (negative) for right (left) polarized waves. 
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Fig. 5. — Exponential of the normalized damping rate of the 
Alfven-whistler at Sj^b = 87° showing cnhacement of the damping 
near the harmonics of ions. 

frequency gaps in the dispersion curves caused by the 
ion Bernstein modes, known to develop in hot plasmas 
at quasi-p e rpren dicular angles of propagation (e.g., Li] 
& Habbal (2001)). The more oblique Bernstein modes 
are the least damped as shown in Fig. [t] (lower panel). 
This result clearly invalidates the magnetosonic disper- 
sion oj ~ kCg of the classical-whistler mode at w > LUci 
found in the previous section from the hot two-fluid the- 
ory. 

Figure [7] shows clearly that, while the Bernstein and 
the fast modes are less damped than the Alfven-whistler 
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Fig. 6. — Dispersion relations of the Alfven-whistler modes (blue) 
and their damping rates (red) for the angles of propagation 80° < 
fkB < 89.99°. 

modes near kpi ^ 1, both the Bernstein and the fast 
magnetosonic modes become heavily damped at kpi > 3. 
The Alfven-whistler modes appear to be the less damped 
ones at scales kpi > 3. Furthermore, among these modes 
the most oblique one is the least damped, which is the 
KAW mode at 6'kB = 89.99° in Fig. [T] 

In Fig. m we tested the validity of these results for dif- 
ferent high f3i values and with = T^. We see clearly 
that the conclusions above remain valid and that the ra- 
tio Ti/Te does not modify much the results but slight 
changes of the dispersion curves near the harmonics of 
ions. This can be seen for instance by comparing Fig. [9] 
(middle panel) to Fig. Is] which was plotted for the same 
$i = 2.9 but with T!j ^oTe. It is important however to 
note that when f3i is becoming smaller (i.e. decreasing 
toward 1) the damping of Alfven-whistler modes by cy- 
clotron resonances become important for some angles of 
propagation, as can be seen in Fig. [9] (top panel). We 
tested other values of /3i (not shown here) and found 
that the conclusions above remain valid. Moreover, we 
found that higher values of f3i allow extending the Alfven- 
whistler modes to even smaller scales (for f3i = (3^ — 25, 
the highest tested value, the modes were found to ex- 
tend up to kpi ^ 30). On the contrary, for (3i 1 
we found that none of the Alfven-whistler modes (when 
80° < 9 < 90°) can propagate at a; > uJd because they 
are strongly damped at a; ^ LUd - 

Now one can ask the question: which of these plasma 
modes, fast, Bernstein, KAW or whistler, is likely to 
carry the energy cascade of turbulence down to the dis- 
sipation scales in the limited range of SW parameters 
studied here? While all these modes might contribute to 
the energy cascade in the SW, based on linear damping 
rates and assuming that (quasi-)linear theory is applica- 
ble to small scale SW fluctuations, we can conclude that 
the KAW branch (w < uJd) is more likely to be observed 
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Fig. 7. — Exponential of the damping rates of the Alfven-whistler, 
fast, and Bernstein modes, normalized to the frequency of each 
wave for angles of propagation 80° < S^b 1^ 89°. 

in the data than the whistler branch (oj > uJd)- How- 
ever, damping rates alone cannot rule out the presence 
of fast or ion Bernstein modes in the transition range 
0.5 < kpi < 3 where they are found to be less damped 
than the Alfven-whistler modes. Why, then, rule out the 
fast and the Bernstein mod es in this ran ge of scales? It 
has been often argued (e.g., Howes [2009 1 that the Bern- 
stein modes are highly electrostatic and therefore cannot 
account for the magnetic spectra observed in the disper- 
sion range in the SW. This is not totally correct as can 
be seen in Fig. [lOj which shows the three components 
of the electric fiela of each mode. With k chosen in the 
XZ plane and considering the high oblique angles stud- 
ied here (80° < 6'kB < 89°) the Ex component (in red 
in the plot) represents essentially the electrostatic field. 
It is instructive to observe that the electromagnetic part 
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Fig. 8. — Dispersion relations of the fast and the ion Bernstein 
modes (blue) and their damping rates (red) for the angles of prop- 
tion 80° < 6»kB < 89°. 



of the Bernstein mode dominates over the electrostatic 
part at scales 0.2 < kpi < 1 and that the electrostatic 
component dominates only at smaller scales. More inter- 
estingly, we observe that the Alfven-whistler modes are 
essentially electrostatic at all scales. This, in fact, is not 
surprising: in the MHD limit the field perturbations of 
the shear Alfven mode SE and S'B and the mean field 
Bq are orthogonal to each other. At quasi-perpendicular 
angles the electric field is therefore necessarily quasi- 
longitudinal (i.e., parallel to k). 

Figure [TO] shows that only the fast mode has a domi- 
nant electromagnetic field up to scales kpi > 3 where it 
becomes comparable to the electrostatic part. The three 
modes have very small parallel electric field (given by 
the Ez component in the plot). Therefore, Fig. 10 shows 
clearly that one cannot rule out any mode only on the ba- 
sis of the strength of its electrostatic component. Instead 
we suggest using the magneti c compressibility (Song et 
ar]|1994 Lacombe et al.||l995 1 , in addition to the damp- 
mg rates discussed above, to rule out the fast and the 
Bernstein modes in SW observations. 

We define the magnetic compressibility by Cb = 
5B\\^ /SB^ (where 5B^ is the total magnetic power). It 
is a well established result that SW magnetic compress- 
ibility is small in the inertial range, i.e., Cb ^0.1, then 
it increases t o Cb ^ 0.4 in the tran s ition and the disper - 
sive ranges ( [Sahraoui et al. 2010a Kiyani et al. 2011). 
In Fig. [iT] we plot the magnetic compressibilities of the 
three modes for the same angles of propagation as before. 
The magnetic compressibility profile of each mode shows 
clearly that the Alfven-whistler modes fit SW observa- 
tions better, at least at the scales kpi > 1, since both 
the Bernstein and the fast modes show the dominance 
of parallel power over perpendicular power (Cb ^ 1)- 
However, it is important to notice in Fig. 11 (top panel) 
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Fig. 9. — Dispersion relations of the Alfven-whistler modes (blue) 
and their damping rates (red) for the angles of propagation 80° < 
^'kB < 89° and for different values of (all have Ti = Te). 
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Fig. 10. — Electric field components (red), Ey (green) and 
Ez (blue) of the Alfven-wfiistler, fast magnetosonic and the ion 
Bernstein modes normalized to their total electric field, computed 
for the angles of propagation 80° < Sj^b < 89°. The E^ component 
(red) is essentially the electrostatic part of the electric field in each 
mode, while E^ (blue) is the parallel component (to Bo) of the 
electric field. 

have the same profile of magnetic compressibility. This 
implies that, at the oblique angles studied here, one can- 
not answer the question as to which branch, KAW or 
the whistler, dominates in the data soley from the mea- 
surement of the magnetic com pressibility. Th i s cont rasts 
with the conclusions given in Gary & Smith (2009). As 
discussed above, the computed damping rates showed 
that the whistler branch is more damped than the KAW 
branch at small scales. 

4. DISCUSSION 

The previous results show that in high f3i plasma and 
at high oblique angles of propagation, there is one solu- 




1 .0 
0.8 
0.6 
0.4 
0.2 
0.0 



T t r I I I r I iiiriiniiii 



Post mognetosonic 



0.01 

1 .0 
0.8 
0.6 
0.4 
0.2 b 
0.0 
0.01 



0.10 1.00 



1 0.00 



Ion Bernstein 




0.10 1.00 
kP: 



1 0.00 



Fig. 11. — Magnetic compressibility of the Alfven-whistler, fast, 
and Bernstein modes for the angles of propagation 80° < ^?kB ^ 
89°. 

tion that extends from the MHD to the electron scales: 
it is the Alfven-whistler mode. The whistler mode dis- 
cussed here develops as a continuation at high frequency 
of the shear Alfven mode and not of the fast magne- 
tosonic mode. This re sult rises some excep tions to the 
general claim made by Howes et al. (2011b) on the ab- 
sence of whistler modes m small scale yw turbulence 
because of the absence of fast magnetosonic modes in 
the inertial range. The fast mode strongly resonates with 
ions near LJci and splits into different ion Bernstein modes 
at higher frequencies. In contrast, the shear Alfven mode 
extends at scales kpi > 1 to frequencies either larger or 
smaller than ujd, depending on the ratio k^\/kj_. The 
same mode can thus be called whistler (w > tOci) or KAW 
{lu < LUci) depending on this ratio. If ku/k± < ^ei then 
the Alfven mode follows the KAW mode with w < uid 
at all scales (even at kpe > 1). For larger values of 
k\\/k±, the Alfven mode develops a whistler branch with 
uj > ujci at a given scale, but the damping rate increases 
as well. Based on the damping rates one would expect to 
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observe more oblique KAW than oblique whistlers. This 
result may be considered as another, rather simpler, al- 
ternative to explain the strong anisotropics {k± » fcy) 
observed in the SW. Indeed, there are several predictions 
from different MHD and kinetic turbulence models, e.g ., 



critical balance conjecture ( Goldreich &: Sridhar ^| 19 95 ) 
that predict stro nger anisotropy at small scales (She- 



balin et al. 19831. Here we argue that the /wear kinetic 
damping oi the Alfven- whistler modes as discussed above 
may as well explain such observed strong anisotropy: the 
more oblique Alfven-whistler modes are indeed the less 
damped ones. 

From this result it appears clear that the question as to 
which branch, KAW or whistler, dominates in the data 
reduces to the question as to how oblique is the wave 
vector k. The answer requires accurate measurement of 
the angle , which in turn requires simultaneous mea- 
surement of ki\ and k±. This is now possible using the 



fc-filtering technique on Cluster spacecr aft data ( Sahraoui 
et al.||2006[|2010a|jNarita et al.||2010a[ ). However, as one 
can see from l'ig.|3j a change in the aiigle of propagation 
from 6'kB = 80° to 0kB ^ 90° will change the physics 
from the whistler branch to the KAW branch. Measur- 
ing the angle 0kB with an accuracy better than 10° is 
quite difficult, even with the powerful multi-spacecraft 
techniques, owing to t he various sources of uncerta. inties 
in the measurements (Sahr aoui et al.|[2003a 2010c). 

When simultaneous measurement of fcy and k± is not 
possible, for instance when only single spacecraf t data are 
available, a nother alternative, evoked briefly in |Sahraoui| 
et al. ( 2009 ), exists to test which mode (whistler of KAW) 
IS present in the data. As we will show below, this 
method is however limited to situations where turbulence 
is not confined to the very oblique angles studied here. 
Let us discuss it in more details here. The general for- 
mula relating the frequency of the wave in the SW rest 
frame w to the measured one onboard the spacecraft uj^c 
is given by 

LOgc = uj + kVsw cos (?kv (14) 

where Vsw is the SW speed forming an angle 0kv with the 
k vector. If the KAW branch (a; < oJd) dominates in the 
dispersive range down to the electron scale, then the Tay- 
lor frozen-in-fl ow assumption should be valid, meaning 
that equation |14| reduces to Usc ^ fcy^^ cos (?kv- Now, 
let us assume furthermore that two breakpoints occur in 
the turbulence energy spectra at the ion and the elec- 
tron gyroscales pi and pe (Sahraoui et al. 2009). As 



shown in Fig. [12] these breakpoints should be observed 
in the spacecraft frame respectively at the Doppler- 
shifted frequencies fp^ = ujp./2Tr = Vsiu cos6'kv/27rpi and 



/p. 



. /27r = y^u) cos0kv/27r/9e. Assuming a similar 



angle 0kv at the ion and the electron scale, the ratio 
between the two frequencies should fulfill the relation 



fh 

Pe 



(15) 



From Fig. [12] one obtains the ratio fp^ / fp. ^ 42 con- 
sidering that it, = Tf.- In the data studied in [Sahraoui] 
et al. ( 2009[ ) equation 15 has been accurately verified: 
with li/l'e ~ 5 the predicted ratio from equation 15 



fp^l fpi ~ 95, which was very close to the ratio be- 
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Fig. 12. — Frequency signatures of the ion and electron gyroscales 
a frozen-in flow approximation (red) and two whistler modes at two 
different angles of propagation are used (T^ = Te). The indicated 
frequencies are given in the text. 



40i?z and ion breakpoint at 



tron breakpoint at j p 
fp^ ^ OAHz. 

If however the relation [l5| is not fulfilled by two ob- 
served breakpoints this would imply that the Taylor as- 
sumption fails, or equivalently, that a dispersive branch 
of the whistler type exists within the data. To show that 
let us now assume that the whistler branch is dominant 
in the dispersive range. In this case the Taylor hypothesis 
will naturally fail (since phase speeds comparable to Vsw 
are present). At the very small scale, i.e., near pe, one can 
even expect that the phase speed of the whistler mode, 
being linear in k (sin ce cu oc fc^), will be larger than Vsw 
as can be seen in Fig.[l2] In this case equation [l4| yields 
Use ^ ^- This means that the relevant temporal scale 
will be given, not by "Doppler-shifting" pe, but by the 
dispersion curve ofthe whistler mode in the plasma rest 
frame. From Fig. 12 one can see that the scale Pe {=dp) 



will yield a different breakpoint frequency depending on 
the angle ^ke^ — 500a;ci and uj2 — lOOOoJci respec- 
tively for 6'kB = 60° and 6»kB = 10°. With /„ - 0.5Hz, 
the breakpoint at the electron scale should thus be ob- 
served onboard the spacecraft at the frequencies 500Hz 
and 250Hz respectively. Thes e frequencies ar e much 



higher than those reported in Sahraoui et al.[ ( |2009[ ) 
which validates the conclusions of that work that tur- 
bulence was more consistent with the KAW modes (or 
quasi-stationary fluctuations) than with whistlers modes. 
At higher oblique angles (6'kB > 60°) one can extrapo- 
late from Fig. [T2[ that the phase speed of the whistler 
branch will become comparable or smaller than Vsw, 
therefore the Doppler-shift cannot be neglected in equa- 
tion [M] An accurate estimation of both terms in the 
right-hand side of equation [T4| is then necessary. As men- 
tioned above, this shows the limit of this simple method 
and emphasizes again the difficulty in distinguishing be- 
tween the two branches, KAW and whistler, of the same 
mode at the very high oblique angles. The estimations 
given here can be tested on SW magnetic power spectra 
that can be measured by the high-time resolution Cluster 
searc h-coils (up to 225H z in the spacecraft frame of refer- 
ence ( Sahraoui et al.|2"00 9, 2010c)) to search for data in- 
tervals when parallel of moderate oblique whistler mode 
turb ulen ce may exist (i.e. intervals not fulfilling equa- 
tion [15|). However, very often, the level of the magnetic 



tween the frequencies corresponding to the observed elec- turbulence in the SW is not sufficiently high and can hit 
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the sensitivity floor of the magnetometers at frequencies 
as low as 20Hz. This hmitation emphasizes the need in 
the future space missions for new search-coils that have 
higher levels of sensitivity in the dispersive and the dis- 
sipation ranges of SW turbulence. 



5. CONCLUSIONS 

Using the hot two-fluid and the Vlasov theories we 
studied key properties of the linear plasma modes un- 
der realistic SW conditions at lAU, namely /3i > /3e ~ 1, 
Ti > Te and focused on high oblique angles of propa- 
gation 89° < ^kB < 90° as observed frequently from 
the Cluster multi-spacecraft data. We discussed the rele- 
vance of each of the main modes (the KAW, the whistler, 
the fast magnetosonic and the ion Bernstein) to carry the 
energy cascade below the ion gyroscale (the slow magne- 
tosonic mode being heavily damped by kinetic effects). 
We showed in particular that the whistler branch devel- 
ops as a continuation at high frequency of the classical 
(shear) Alfven wave known in the MHD limit, and not 
of the fast magnetosonic mode. The Alfven mode ex- 
tends indeed at scales kpi ^ 1 to frequencies either larger 
or smaller than ujd, depending on the anisotropy fcy/fcj^. 
The same mode is thus called whistler (w > LOd) or KAW 
(w < LOci) depending on the anisotropy ratio. 

Unlike the two-fluid model, the Vlasov theory shows 
that the fast magnetosonic mode undergoes strong cy- 
clotron damping at a; ~ ujd and splits up into ion Bern- 
stein modes at a; > Wci, in agreement with the analysis 
of the wave polarization. We showed also that ruling 



out the fast and the Bernstein modes from the cascade 
process in the dispersive range based on their electro- 
static nature is incorrect, since the Alfven-whistler modes 
are shown to be electrostatic as well at high oblique an- 
gles. We suggest rather that one must combine mag- 
netic compressibility and damping rates to rule out those 
modes: both were shown to be highly compressible (i.e., 
By >> B']_) in disagreement with typical SW observa- 
tions of magnetic field data. In addition, those modes 
were shown to be heavily damped in the dispersive range 
as compared with the Alfven-whistler mode. We showed 
finally that the modes on the whistler branch are more 
damped than those of the KAW branch, which makes 
the latter a more relevant candidate to carry the energy 
cascade down to electron scales. We also demonstrate a 
way to test the presence of either branch in SW space- 
craft data and discussed the limitation of the existing 
spacecraft observations that need to be kept in mind in 
the planning of the future space missions. 

Finally, it is important to keep in mind that the con- 
clusions of this work may not be applicable to all SW ob- 
servations because the present study is fully linear (one 
can question the relevanc e of linear theor i es to S W tur- 
b ulence , see for instance Dmitruk et al. (20091; Howes 
|](2009)) and restricted to the range ot parameters dis- 
cussed above. Moreover, the role that can be played 
by several plasma instabilities (both at ion an d electron 
scales ) as re p orted in several s pace observations Sahrao ui 



et al.| ( |2006[ ); rHellinger et al.| ( |2006[ ); |Bale et al. 

not addressed here. 
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